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BR1 



DB 540 9 ff / 210502 



™ BE 2? 6 3UIN 2003 
UEU 75 IN PI PARIS B 

0307848 



N° D'ENREGISTREMENT 
NATIONAL ATTRtBUE PAR LINPI 

I DATE DE DEPOT ATTR1BUEE 
I PAR L'INPI 



2 7 JUIM 2003 



I Vos references pour ce dossier 

\(facultaiif) B 143643 JCI DP 2528 



Confirmation d'un depot par telecopie 



0 NATURE DE LA — ^ 



Demande de brevet 



Demande de certificat d'utilite 



J NOW! ET ADRESSE DU DEMANDEUR OU DU MANDATAIRE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 

■ 

BREVATOME 

3, rue du Docteur Lancereaux 
75008 PARIS 
422-5 S/002 



PI N° attribue par I'lNPI a la telecopie 



Cochez Tune des 4 cases suivantes 



El 



□ 



Demande divisionnaire 

Demande de brevet initiate 
ou detnatide de certificat dhdilM initiate 



Transformation d'une demande de 
brevet europeen Demande de brevet initiate | N 



□ 

N° 
□ 



Date 
Date 



1 I I I I ■ ' I I 



Date I i 1 , I 



l I ' 



I jgcj TITRE DE L' INVENTION (200 caracteres ou espaces maximum) 



j Q DECLARATION DE PRIORITY 
OU REQUETE DU BENEFICE DE 
LA DATE DE d£POT D'UNE 
DEMANDE ANTERIEURE FRANQA1SE 



Pays ou organisation 
Date hhl' 



Pays ou organisation 
Date | t I i I I i 



Pays ou organisation 
Date 11)1 



□ S'il y a d'autres priorites, cochez la case et utilisez I' imp rim 6 «Suiten 



DEMANDEUR (Cochez Tune des 2 cases) 



[71 r^6te%a1e O Person ne physique 



Norn 

ou denomination sociale 



COMMISSARIAT A L'ENERGIE . ATOMIQUE 



Prenoms 



Forme juridique 



Etablissement Public de Caractere Scientifique, Technique et Industriel 



N° SIREN _ 
Code APE-NAF 



I i i i i i i » i 1 



j L 



J 



Domicile 
ou 



Rue 



3 1-33, rue de la Federation 



Code postal et ville 



Pays 



Nationality 



N° de telephone (facultalif) 



Adresse electronique (facidiatif) 



|7 5|7:5,2|PARIS 15eme 
FRANCE 



francaise 



N° de telecopie (facullatif) 



\E\ S'il y a plus d'un demandeur, cochez la case et utilisez Pirn prime «Sulte» 



Remplir imperativement la ^_page 



INPI 



■ IKSTITUI 
NATIONAL OE 
LA PROPB1ETE 



INOOSTRIBLLI 

26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 33 (1) 53 04 53 04 Telecopie : 33 (1) 42 94 86 54 



REMISE OES PIECES 
DATE 



jReserveat'lNPl} 



t LMNPl 

DATE DE DIiFOT ATTRIBUTE % \\ <S 



UEU 

N° D'ENREGJSTREMENT 
NATIONAL ATTR)BU£ PAR L'fNPl 



PAR L'tNPf 
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(facultatij) B 14364.3/JCI DD 2528 
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get imprime est a remplir lisiblement a I'encre noire 




N° 11354*03 



DB 540 98 / 210502 



Q NOR/1 ET ADRESSE DU DEMANDEUR OU DU MANDATAtRE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSfc 

BREVATOME 

3, rue du Docteur Lancereaux 
75008 PARIS 
422-5 S/002 



Confirmation d'un depot par telecopie 



□ N° attribue par I'lNPI a la telecopie 



Demande de brevet 



® 



Demande de certificat d'utilite 



□ 



Demande division naire 

Demande de brevet initiate 
ou demande de certificat d'idilit& initiate 



□ 

N° 



Date Li I i 1 i I i I 

Date 1 i 1 i I i i i 1 



Transformation d'une demande de 
brevet europeen Demande de brevet initiate 



□ 
N° 



■■1.1 



1 i i 



El T1TRE DE L'INVENTION (200 caracteres ou espaces maximum) 

PROCEDE DE RECONSTRUCTION D'UNE IMAGE TOMOGRAPHIQUE PAR UNE METHODE 
ANALYTIQUE COMPRENANT UNE MODELISATION AMELIOREE DU MOUVEMENT DE 
L'OBJET. 



DECLARATION DE PRIORITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANTERIEURE FRANQAISE 



Pays ou organisation 

i i i i i i i i 



Pays ou organisation 

Date I | | | | i | | | N° 
Pays ou organisation 

Date I i 1 i I i i i 1 N° 

□ SMI y a gauges priorftes, cochez la case et utilisez I'imprime «Suite» 



Norn 

ou denomination sociale 



COMMISSARIAT A V ENERGIE ATOMIQUE 



Prenoms 



Forme juridique 



Etablissement de caractere Scientifique, Technique et Industriel 



N° SIREN 



Code APE-NAF 



1 1 I I I M I I I 



I I I I I 



Domicile 
ou 



Rue 



31-33 rae de la Federation 



Code postal et ville 



Pays 



17 15,7.5,2 1 PARIS 15eme 



FRANCE 



Nationality 



FRANCAISE 



N° de telephone (faadtatif) 



N° de telecopie (faadtaHf) 0 



Adresse electronique (facutlatif) 



I I S'il v a nluft rf'nn ffamanHAiir. r.netha-r la csca ot it+Sli<so7 Vimnrlma o.Qu(tAn 




NATIONAL Ot 



BREVET D'INVENfl 
CERTIFICAT D'UTHJl 



^1 
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BR2 



KfSrSluiN 2003 

ueo 75 IN PI PARIS B 

030784S 



IrteservSa riNPIf 



N« D'ENREGISTREMENT 
NATIONAL ATTRIBUE PAR L'INPI 



DB 540 W/ 210502 







Norn 


LEHU 


Prenom 


Jean 


Cabinet ou Societe 


BREVATOME 


N °de pouvoir permanent eVou 
de lien contractuel 


PG 7068 


Adresse 


Rue 


j 5 rue qu JL/ocicur j-fduocrcaUA 


Code postal et ville 


17 i5 iO iO 18 J PARIS 


Pays 


FRANCE j 


N° de telephone (factdtoti/) 


01 53 83 94 00 


N° de telecopie (faadtattj) 


01 45 63 83 33 


Adresse Electronique (facuiiattj) 


brevets.patents@brevalex.com 


m mffiteuKW" ' ; ^^^ ^ 


Les fhven&uresorSf n^essalreme^ itie's personnes pKysiqties : 7 ' ,; "f : " : * 


Les demandeurs et les inventeurs 
sont les memes person nes 


□ Oui 

fXl Non : Dans ce cas remplir le formulaire de Designation d'inventeur(s) ? : 


Wt raSppoitcW^ 


Uniquement £our iihe demanoe de brevet (y 


Etablissement tmmediat 
ou etablissement differe 


a 
n 


Paiement echelonne de la redevance 

{en deux versements) 


Uniquement pour les personnes physiques effectuant elles-memes leur propre depot 

UOui 

□ Non 


O REDUCTION DU TAUX 
DES REDEVANCES 


Uniquement pour les personnes physiques 

I | Requise pour la premiere fois pour cette invention (/oindtv tin avisde non-imposition) 
1 I Obtenue anttrieurement a ce depot pour cette invention (joindre une cop/e de /a 
d€a'siond I ad7nissional , assistencegratuiteouit?diguersartf 1 > i i i 1 


SEQUENCES DE NUCLEOTIDES 
ET/OU D'ACIDES AMINES 


1 1 Cochez la case si la description contient une liste de sequences 


Le support electronique de donnees est joint 

La declaration de conformite de la llste de 
sequences sur support papier avec le 
support electronique de donnees est jointe 


□ □ 


Si vous avez utilise* I'imprime «Suite», 
indiquez le nombre de pages jointes 


i 



SIGNATURE DU DEMANDEUR 

OU DU W1ANDATAIRE 

(Noni et qualite du signataire) 



J. LEHU 




OU DE L'INPI 




La loi n°78-17 du 6 Janvier 1978 relative a rinformatique, aux fichiers et aux libertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de riNPI. 



regue te ^3/u//u3 




INST HOT 
NATIONAL DE 
LA FROPRIETE 
IMO US TBI ELL G 



26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 33 (1) 53 04 53 04 Telecople : 33 (1} 42 94 86 54 



REMISE DES PIECES 
DATE 

UEU 

N° D'ENREGISTREMENT 
NATIONAL ATTR1BUE PAR L'INPI 



- I Reserve a llNPil 



B&EVET INVENTION 
CERTIFICAT D'UTILITE 

Code de la propriety intellectueile - Uvre VI 

REQUETE EN DELIVRANCE 

Page suite N° .../.. 




N° 11354*03 




Cet imprime est a remplir lisiblement a 1'encre noire 



DB 829 ©W/ 010702 



Vos references pour ce dossier {faculiaiij] 


B 143643/JCI DD 2528 


Q DECLARATION DE PRIORITE 
OU REQUITE du BENEFICE DE 
LA DATE DE DEPOT D'UNE 
DEMANDE ANTERIEURE FRANQAISE 


Pays ou organisation 

n*te 1 i 1 i 1 i i i 1 n° 

Pays ou organisation 

iii i 

natP 1 1 1 1 1 1 1 1 1 N° 

Pays ou organisation 

rtata \ I | 1 1 1 1 j 1 MO 

uaie ii i i i I i i i in 






Nom 

ou denomination sociaie 


CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 


Prenoms 




Forme juridique 




N° SIREN 


1 i i i i i i i i 1 


Code APE-NAF 


l.i.l 


Domicile 
ou 

siege 


Rue 


3 rue Michel Ange 


Code postal et vllle 


I7i5i7i9i4I PARIS CEDEX 16 


Pays 




Nationality 




N° de telephone {facultatifl 




N° de telecopie Ifacultatift 




Adresse electronique {faculiatift 








Nom 

ou denomination sociaie 




Prenoms 




Forme juridique 




N° SIREN 




Code APE-NAF 


i.i.i 


Domicile 
ou 

siege 


Rue 




Code postal et vllle 


i * i i i i 


Pays 




Nationality 




N° de telephone tfacultatifi 




N° de telecopie \facultalif\ 




Adresse electronique ifaciiltaHj] 




2J SIGNATURE DU DEMANDEUR 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 

J.LEHU 


VISA DE LA PREFECTURE 
OU DE L'INPI 



La ioi n°7S-17 du 6 janvier 1978 relative a I'lnformatique, aux fichiers et aux Iibert6s s'applique aux rSponses faites a ce formulafre. 
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i net i tut 
national oh 
la propriete 
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26 bis, rue de Saint P&ersbourg 
75800 Paris Cedex 08 

Telephone : 33 (1) 53 04 53 04 Telecopie : 33 {I) 42 94 86 54 



BREVET D'U^/EIMTION Kwpi 
CERTIFICAT IH'ILITE 

Code de la proprtete intellectuelle - Livre VI 

REQUETE EN DELIVRANCE 



SSfSWUIN 2003 
ueu 75 INPI PARIS B 

0307848 



) Reserve a TINPlt 



Page suite N° \ . . f\ 



BR 



,/SUlTE 



N° D'ENREGISTREMENT 
NATIONAL ATTRIBUE PAR UINPI 



Cet imprime est a remplir lisibfement a I'encre noire ob 829 © w / 010702 



Vos references pour ce dossier Ifacultatij) 




□ DECLARATION DE PRIORITE 
OU REQUETE DU BENEFICE DE 
LA DATE DE D£p6T D'UNE 
DEMANDE ANTERIEURE FRANpAlSE 


Pays ou organisation 

nate 1 1 1 - 1 1 1 1 1 N° 
Pays ou organisation 

n a to 1 1 1 1 1 1 l 1 1 N° 
Pays ou organisation 

Date 1 L 1 N° 




Q0 Personhe mo^le; : * ' : v J' ,.E3-?J?!S?K!!?.?1 


Norn 

ou denomination sociale 


CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 


Prenoms 




Forme juridique 




N° SIREN 


1 1 1 1 i 1 1 1 1 '■ 


Code APE-NAF 




Domicile 
ou 

siege 


Rue - 


3, rae Michel Ange 


Code postal et ville 


17 -5 7 s 9 41 PARIS CEDEX 16 C 


Pays 


FRANCE 


Nationality 


francjaise 


N° de telephone {facutiatifl 


r 


N° de telecopie {facultatif] 




Adresse electronique ifacultatij) 






□ Personiie.fflf&e* " * r Q Perspftne physique 


Norn 

ou denomination sociale 




Prenoms 




Forme juridique 




N° SIREN 


l 1 1 1 1 1 11! : 


I Code APE-NAF 


1 1 1 1 1 


Domicile 
ou 

siege 


Rue 




Code postal et ville 




Pays 




Nationality 




N° de telephone \facultalij\ 




N° de telecopie {facullatifi 




Adresse electronique \facultatifi 




M SIGNATURE DU DEMANDEUR > 
OU DU MANDATAIRE / «d 

(Norn et qualite du signataire) 1 \ s 

J. LEHU lA^/ ^ 


VjrtSA DE LA PRipECTURE 



La loi n°78-17 du 6 janvier 1978 relative a I'informatique, aux fichiers et aux liberies s'applique aux reponses faltes a ce formulaire. 
EHe garantit un droit d'acces et de rectification pour les donriees vous concernant aupres de I'INPI 
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CERTIFICAT D'UTILITE 



REQUETE EN DELIVRANCE 
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m 



REMISE OES PIECES 
DATE 

UEU 

N* D'ENREGISTREMENT 
NATIONAL ATTRIBUE PAR LINPI 



\ Reserve a 11NPI h 







Nom 


J^JbxiU 


Prenom 


Jean 


Cabinet ou Soctete 


BREVATOME 
422.5/S002 


N °de pouvoir 
de lien contra< 


permanent et/ou 
;tuel 


7068 du 12.06.98 


Adresse 


Rue 


3, rue du Docteur Lancereaux 


Code postal et ville 


17 i5 iO iO 18IPAPTS 


Pays 


FRANCE 


N° de telephone Cfacultattf) 


01 53 8*3 04 on 

VJA JJ OJ y-r \J\J 


N° de teiecopie <faculfat0 


01 45 63 83 33 


Adresse electronique (fauifaftfl 


1"*"rfiVtf*t<3 llfltPTltQf^hrAvx/alAY r*rvm 


Les demandeurs et les inventeurs 
sont les memes personnes 


Clf Oui ^ 

13 Non : Dans ce cas remplir le formulaire de Designation d'inventeur(s) 


Etablissement immediat 
ou 6tablrssement dlffere 


□ 


Paiement echeionne de ia redevance 

(en deuxvetsetneitte) 


Uniquement pour les personnes physiques effectual elles-memes leur propre depot 

LJ Oui 
□ Non 


B REDUCTION DU TAUX 
DES REDEVANCES 


Uniquement pour !es personnes physiques 

□ Requlse pour la premiere fois pour cette Invention (/otrufrvrmaiwdemmmjHmWoxJ 

□ Obtenue anterieurement a ce ctepot pour cette invention tfomtfre une copie de /a 
Mcfsiond'admissiond/'twisfaftcegt-afyfe 1 i | | | | 


EE SEQUENCES DE NUCLEOTIDES 
ET/OU D'ACIDES AMINES 


O Cochez la case si la description content une liste de sequences 


Le support electronique de donnees est joint 

La declaration de conformite de la liste de 
sequences sur support papier avec le 
support electronique de donnees est jointe 




S! vous avez utilise Fimprime «5uite», 
indiquez le nombre de pages Jointes 




EH SIGNATURE DU DEMANDEUR 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 

J. LEHU 


VISA DE LA PREFECTURE 
OU DE L'INPI 



La loi n°78.17 du 6 janvier 1978 relative a I'informatique, aux fichiers et aux liberies s r applique aux reponses faites a ce formulaire. 

EMe ffarantlt lin rlmft H'arrA*: a* Ha rArtffi/»a*iAn nmir lac Hnnnoot- wnur ^^r,^ rt ^„^„+ „.,^x^ ninint 



PROCEDE DE RECONSTRUCTION D 1 UNE IMAGE TOMOGRAPHIQUE PAR 
UNE METHODE ANALYTIQUE COMP RENANT UNE MODELI SATION 
AMELIOREE DU MOUVEMENT DE L 1 OB JET • 

5 DESCRIPTION 

Le sujet de cette invention est un proced€ 
de reconstruction d'une image tomographique, du genre 
analytique et dans lequel on recourt §l une modelisation 

10 perfect ionnee du mouvement de l'objet afin de r^duire 
les artefacts de l 1 image. 

. La reconstruction d 1 images en tomographie 
implique 1 'utilisation d'un rayonnement traversant 
l'objet. Dans de norabreux cas, le rayonnement est un 

15 rayonnement d f irradiation partiellement attenue par les * 
points de I'objet qu'il traverse / il arrive aussi quef; 
les lignes de rayonnement soient fictives et 

correspondent simplement k des lignes de collimation k.\ 

-■I 

des detecteurs, qui enregistrent une emission de V 
20 particules dans I'objet. Les premiers de ces procedes 
sont les procSdes radiographiques , et les autres des 
proced6s Smissifs, reserves plut6t aux etres vivants et 
dans lesquels le rayonnement est produit par un corps 
emissif qu'ils ont absorbs auparavant . Quoique les deux 
25 families de procedes soient completement differentes, 
la plupart des procedes de reconstruction d' images 
s'appliquent aux deux, et c'est le cas aussi de 
1 1 invention. 

Un rayonnement divergeant est souvent 
30 prefers pour englober plus facilement l'objet £tudi<§. 
Dans le cas d ! un rayonnement d f irradiation attenue par 
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reyue le <^/u//ua 



PROCEDE DE RECONSTRUCTION D'UNE IMAGE 
TOMOGRAPHIQUE PAR UNE METHODS ANAL YT I QUE COMPRENANT UNE 
MODELI SATION AMELIOREE DU MOUVEMENT DE Ii 1 OBJET • 

5 DESCRIPTION 



Le sujet de cette invention est un proced6 de 
reconstruction d'une image tomographique , du genre 
analytique et dans lequel on recourt a une mod^lisation 

10 perf ectionnee du mouvement de I'objet afin de reduire 
les artefacts de 1 1 image . 

La reconstruction d ! images en tomographie 
implique 1 1 utilisation d 1 un rayonnement traversant 
I'objet. Dans de notnbreux cas, le rayonnement est un 

15 rayonnement d f irradiation partiellement attend par les 
points de I'objet qu f il traverse ; il arrive aussi que 
les lignes de rayonnement soient fictives et 
correspondent simplement a des lignes de collimation 
des detecteurs, qui enregistrent une emission de 

20 particules dans I'objet. Les premiers de ces procedes 
sont les procedes radiographiques, et les autres des 
procedes £missifs , , reserves plutot aux etres vivants et 
dans lesquels le rayonnement est produit par un corps 
£missif qu'ils ont absorbe auparavant . Quoique les deux 

25 families de procedes soient complStement diffSrentes, 
la plupart des procedes de reconstruction d f images 
s'appliquent aux deux, et c'est le cas aussi de 
1 ' invention. 

Un rayonnement divergeant est souvent prefere 
30 pour englober plus facilement I'objet etudie. Dans le 
cas d'un rayonnement d' irradiation att^nu^ par I'objet, 



• 2 • 



1' objet, on y parvient en disposant une source emissive 
a un foyer ponctuel pres de I 1 objet et un reseau de 
detecteurs du c6te oppose, de l 1 objet, qui sont tous 
collimatSs vers la source ; dans le cas d'un 
5 rayonnement 6mis par I 1 objet, un rgseau similaire de 
detecteurs est mis en collimation vers un foyer qui 
correspond g^ometriquement a la source ponctuel le 
prec£dente, mais h aucun objet materiel. Dans tous les 
cas, le foyer et le reseau de d£tecteurJsont deplac<§s J*"* 

10 autour de 1* objet en en prenant des vues successives, 
et pour chacune des vues, les detecteurs mesurent des 
sommes de la propriety d 1 attenuation ou demission le 
long de lignes de collimation, qu ! on appelle des 
projections de 1 ' image de 1' objet. Quand un notnbre 

15 suffisant de vues a ete pris, on dispose d f un grand 

•V. 

notnbre de projections s 1 entrecroisant a travers " 
1» objet. Le probldme mathematique dit d 1 inversion du S 
systeme eu de — -retro-projection permet d'obtenir la 
propriete d 1 attenuation ou d 1 emission, qui sert a 

20 former l 1 image, a chacun des points de l f objet S partir ? 
des sommes de proprietes sur les projections. II peut 
etre represents par tin systdme lin^aire d 1 equations ou 
les valeurs connues au depart sont les mesures prises 
par les detecteurs aux differents emplacements et les 

25 inconnues sont les valeurs de la propriete " aux 
differents points de l 1 objet. Certaines methodes de 
reconstruction dites algebriques reviennent 

ef f ectivement & une inversion de ce systeme ; il existe 
cependant d'autres methodes dites analytiques, ou la 

30 valeur de la propriete a chaque point est calculee 

I 
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on y parvient en disposant une source emissive a. un 
foyer ponctuel pres de l"objet et un reseau de 
detecteurs du c6te oppose de l'objet, qui sont tous 
collimates vers la source ; dans le cas d'un 
rayonnement emis par l'objet, un reseau similaire de 
detecteurs est mis en collimation vers un foyer qui 
correspond geometriquement a la source ponctuel le 
precedente, mais a aucun objet materiel. Dans tous les 
cas, le foyer et le reseau de detecteurs sont deplaces 
autour de l'objet en en prenant des vues successives, 
et pour chacune des vues, les detecteurs mesurent des 
sommes de la propriety d' attenuation ou d' emission le 
long de lignes de collimation, qu'on appelle des 
projections de 1 ' image de l'objet. Quand un nombre 
suffisant de vues a ete pris, on dispose d'un grand 
nombre de projections s ' entrecroisant a travers 
l'objet. Le probleme mathematique dit d' inversion du 
systeme permet d'obtenir la propriety d' attenuation ou 
d' emission, qui sert a former 1' image, a chacun des 
points de l'objet a partir des sommes de proprietes sur 
les projections. II peut etre represente par un systeme 
lineaire d' equations ou les valeurs connues au depart 
sont les mesures prises par les detecteurs aux 
differents emplacements et les inconnues sont les 
valeurs de la propriety aux diff6rents points de 
l'objet. Certaines methodes de reconstruction dites 
algebriques reviennent ef f ectivement a une inversion de 
ce systeme ; il existe cependant d'autres methodes 
dites analytigues, ou la valeur de la propriety a 
chaque point est calculee directement a partir d'une 
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directement a partir d'une corabinaison math^matique des 
project ions . L 1 invention est I'une d'entre elles . 

Un domaine privilegie de la tomographie est 
l 1 etude des Stre vivants pour la medecine. Corame 

5 plusieurs vues successives doivent §tre pr^ses^l f ( ^ e1 ^-, n 

£tudi£ est susceptible de se deplacer entre elles. Les A"**&*J$^ 
diff<Srentes vues repr^sentent done des £tats differents 
de l'objet, et des artefacts de reconstruction ne 
peuvent qu 1 apparaitre sur 1 1 image reconstruite . 

10 Des procedis pour Schapper a ce probleme 

permanent consistent & prendre toutes les vues a la 
fois avec autant de sources et de reseau de detecteurs, 
ou & ne prendre des vues de I'objet qu'a des etats * 
identiques de celui-ci, ce qui est possible si son 

15 mouvement est periodique comme le sont certaines ■* 
activites physiologiques comme les battements du cceur -tf 
ou la respiration ; le premier de ces proc^des est 
cependant dispendieux et le second est d^licat & mettre 
en oeuvre convenablement. 

20 C'est pourquoi on a aussi entrepris 

d'incorporer des modeles numeriques du mouvement de 
I'objet^ans le proced£ de reconstruction, Le raodele 
revient g6neralement a une collection de champs de 
deplacements des points de l'objet aux instants 

25 respectifs ou les vues sont prises. II peut §tre <§tabli 
par d'autres mesures ou par certaines hypotheses. 
Plusieurs genres de modeles ont dejS. ete proposes, mais 
ils presentent en general des insuf f isances . C'est 
ainsi que le brevet US A 5 287 276/ applique a I'gtude 

3 0 de la cage thoracique, considere les translations et 
les dilatations de son contenu, mais ce modele est 
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combinaison mathematique cies projections, L f invention 
est l f une d'entre elles. 

Un domaine privil<§gi€ de la tomographie est 
l'€tude des etre vivants pour la medecine. Comma 
5 plusieurs vues successives doivent §tre prises, l'objet 
etudiS est susceptible de se deplacer entre elles selon 
une loi generale devolution. Les differentes vues 
representent done des gtats diff brents de l'objet, et 
des artefacts de reconstruction ne peuvent 

10 qu 1 apparaitre sur 1 1 image reconstruite . 

Des procedes pour echapper & ce problSme 
permanent consistent a prendre toutes les vues k la 
fois avec autant de sources et de reseau de detecteurs, 
ou a ne prendre des vues de l'objet qu f & des etats 

15 identiques de celui-ci, ce qui est possible si son 
mouvement est pgriodique • comme le sont certaines 
activites physiologiques comme les battements du cceur 
ou la respiration ; le premier de ces procedes est 
cependant dispendieux et le second est d^licat a mettre 

20 en oeuvre convenablement . 

C'est pourquoi on a aussi entrepris 
d'incorporer des modeles numeriques du mouvement de 
l'objet refl€tant la loi generale d 1 evolution de 
l'objet dans le proc£d£ de reconstruction. Le modele 

25 revient generalement a une collection de champs de 
deplacements des points de l ! objet aux instants 
respectifs ou les vues sont prises, II peut etre £tabli 
par d'autres mesures ou par certaines hypotheses. 
Plusieurs genres de modules ont dejd. gte proposes, mais 

30 ils presentent en general des insuff isances. C'est 
ainsi que le brevet US A 5 287 276, appliqu6 a 1' etude 



insuffisant pour d'autres mouvements comme celui du 
coeur, qui comprend une torsion. De plus, les jeux de 
projections utilises pour appliquer la formule 
d' inversion ne sont pas connus rigoureusement, ce qui 
5 peut permettre de nouveaux artefacts. 



avantages d'utiliser un modele plus dev-e-lopp6 mais 
aussi tres simple du mouvement de l T objet, comprenant 
de nouvelles classes de mouvements, et de recourir a 
10 une formule d' inversion nouvelle et parfaitement 
exacte . 



de projections necessaires et seulement necessaires a 
la reconstruction des points de 1' image . 



- 1 1 explication de 1 1 invention : 

- la figure 1 illustre la prise des mesures 
de l'objet etudie, 

- la figure 2 illustre une representation 
25 equivalente l'objet €tant maintenu & un etat de 

reference , 



Le procede de 1 1 invention 




les 



De plus, il est facile d'obtenir les jeux 




- les figures 3a, 3b, et 3c illustrent un 
cas complexe de deformation de l'objet, 

- la figure 4 analogue a la figure 2 
3 0 dans le cas complexe, 
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de la cage thoracique, considdre les translations et 
les dilatations de son contenu, mais ce modele est 
insuffisant pour d ! autres mouvements comme celui du 
coeur, qui comprend une torsion. De plus, les jeux de 
5 projections utilises pour appliquer la formule 
d 1 inversion ne sont pas connus rigoureusement, ce qui 
peut permettre de nouveaux artefacts. 

Le procede de 1 1 invention pr6sente les 
avantages d'utiliser un modele plus €labor£ mais aussi 
10 tr^s simple du mouvement de l'objet, comprenant de 
nouvelles classes de mouvements, et de recourir S. une 
formule d 1 inversion nouvelle et parfaitement exacte. 

De plus, il est facile d'obtenir les jeux de 
projections necessaires et seulement ngcessaires k la 
15 reconstruction des points de l 1 image. 

Des generalisations a des situations plus 
complexes sont aussi proposees . 

Les figures suivantes sont introduites pour 
l 1 explication de l 1 invention : 
20 - la figure 1 illustre la prise des mesures de 

l 1 ob jet §tudie, 

- la figure 2 illustre une representation 
gquivalente l'objet Stant maintenu a un etat de 
reference, 

25 - les figures 3a, 3b, et 3c illustrent un cas 

complexe de deformation de l'objet, 

- la figure 4 est analogue a la figure 2 dans 
le cas complexe, 

- la figure 5 illustre une prise de mesures 
3 0 avec un rayonnement conique, 



- la figure 5 illustre une prise de 



mesures avec un rayonnement conique^ 

On aborde maintenant les figures . La 
premiere d'entre elles illustre la situation classique 
en tomographie d'une source S de rayonnement et d'un 
r£seau de d^tecteurs D mobile de concert sur des 
trajectoires circulaires concentriques k des positions 
opposees autour d T un objet E a etudier . Des rayons R 
joignent la source S aux detecteurs respect if s du 
r6seau D . On considere ici un probleme geometrique 
plan, qui correspond aux conditions de mesures dites en 
€ventail parallSle o\x I 1 etude de 1' objet E est faite 
par tranches superpos^es. Les procedes impliquant un 
rayonnement conique sont aussi usuels . et seront 
examines plus loin, L r invention a ete congue pour des 
raypnnements divergent s en general , qui englobent- aussi 
ces derniers . 

La trajectoire T de la source S et 
1 1 ouverture du faisceau de rayonnement sont 
g^neralement choisies pour que celui-ci englobe toute 
la section de l f objet E. 

Considerons un point P de 1' objet E. Un 
rayon unique provenant de la source S le traverse quand 
une vue est prise, et d'autres rayons R2, R3, etc* le 
traversent de meme pour d'autres positions S2, S3, etc, 
de la source S quand d'autres vues sont prises. Ces 
rayons passant par le point P sont utilises pour 
determiner 1 1 image du point P dans le proced£ 
d 1 inversion. II est. utile que les rayons consideres 

soient orientes dans des directions aussi variees que 

i 
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- et les figures 6, 7 et 8 sont des 
organigrammes de trois modes de realisation de 
1 1 invention. 

On aborde maintenant les figures. La premidre 
5 d'entre elles illustre la situation classique en 
tomographie d'une source S de rayonnement et d'un 
reseau de detecteurs D mobile de concert sur des 
trajectoires circulaires concentriques a des positions 
opposees autour d'un objet E £ etudier. Des rayons R 

10 joignent la source S aux detecteurs respectifs du 
reseau D. On considdre ici un probldme geomStrique 
plan, qui correspond aux conditions de mesures dites en 
eventail parallele oti 1» etude de 1' objet E est faite 
par tranches superposees. Les procedes impliquant un 

15 rayonnement conique sont aussi usuels et seront 
examines plus loin. L 1 invention a ete congue pour des 
rayonnement s divergent s en general, qui englobent aussi 
ces derniers. 

La trajectoire T de la source S et l'ouverture 

20 du faisceau de rayonnement sont gen^ralement choisies 
pour que celui-ci englobe toute la section de 1" objet 
E. 

Consid<§rons un point P de 1 ■ objet E. Un rayon 
unique provenant de la source S le traverse quand une 

25 vue est prise, et d'autres rayons R2, R3, etc. le 
traversent de m§me pour d'autres positions S2, S3, etc. 
de la source S quand d 1 autres vues sont prises . Ces 
rayons passant par le point P sont utilises pour 
determiner 1 1 image du point P dans le procede 

30 d' inversion. II est utile que les rayons consideres 
soient orientes dans des directions aussi variees que 



6 

possible,. ^-esJs^^^Ee-jq^^ sot iree S accomplisse au- 

Considerons maintenant que l'objet E se ^ 

deplace et se dSforme. Le point P se deplace a 

5 l f int£rieur de la trajectoire T et les rayons R2, R3, 

etc • a employer passent par des positions respectives 

P2, P3, etc* distinctes du point P d T origine. Le 

probl^me de reconstruction a r^soudre est identique a 

la configuration fictive de la figure 2, o\i le point P 

10 a- ete pris comme reference immobile et ou tout se passe 

comme si la source S suivait une trajectoire T 1 de 

forme irr^guliere. On verra que cette representation 

fictive du probldme aide §l sa comprehension. D'apres 

l 1 invention, le moddle de deplacement et de deformation 

15 qu'on choisit pour l'objet E est af f ine ,^ conf orm^ment §L 

I 

la formule (1) r ~ ^ A o V ^ ^ o - 

20 oil x 0 est la position vectorielle du point P 

par rapport & une reference telle que le point O, §l 
1 1 instant de reference choisi pour mener les calculs et 

la reconstruction de l 1 image, x est la position du 
point P a un autre instant et notamment un instant de 
2 5 prise de vue, et les coefficients a et b dependent du 
temps. X est un parametre general du procSde, qui varie 
entre une valeur minimale X min et une valeur maxima le 
X max ; la position de la source S est reperee par le 
vecteur a. {X) et le temps par une fonction monotone de 
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possible , de fagon S decrire toutes les directions de 
droite possibles pour chaque point de l f objet E. 

Considerons maintenant que l'objet E se deplace 
et se deforme. Le point P se deplace k 1 ■ intSrieur de 
5 la trajectoire T et les rayons R2, R3 7 etc. & employer 
passent par des positions respectives P2, P3, etc. 
distinctes du point P d'origine. Le probleme de 
reconstruction a rSsoudre est identique a la 
configuration fictive de la figure 2, ou le point P a 

10 et6 pris comme reference immobile et ou tout se passe 
comme si la source S suivait une trajectoire T 1 de 
forme irr^guliere. On verra que cette representation 
fictive du probleme aide k sa comprehension. D'apres 
l 1 invention, le moddle de deplacement et de deformation 

15 qu'on choisit pour l'objet E est affine, composee par 
exemple d'une combinaison variable en cours 
d 1 acquisition de translations , de rotations et 
d f homoth<§ties, conformeraent a la formule (1) 

oH x 0 est la position vectorielle du point P par 
rapport a une reference telle que le point O, k 
I 1 instant de reference choisi pour mener les calculs et 
25 la reconstruction de l 1 image, x est la position du 
point P I un autre instant et notamment un instant de 
prise de vue, et les coefficients a et b dependent du 
temps. X est un paramStre gSn^ral du precede, qui varie 
entre une valeur minima le X min et une valeur maximale 



6 



possible, de fagon a decrire toutes les directions de 
droite possibles pour chaque point de l'objet E. 

Considerons maintenant que l f objet E se deplace 
et se deforme. Le point P se d6place a 1 f interieur de 
la trajectoire T et les rayons R2, R3, etc. a employer 
passent par des positions respectives P2, P3, etc. 
distinctes du point P d'origine. Le probleme de 
reconstruction a resoudre est identique a la 
configuration fictive de la figure 2, ou le point P a 
ete pris comme reference immobile et ou tout se passe 
comme si la source S suivait une trajectoire T 1 de 
forme irreguliere. On verra que cette representation 
fictive du probleme aide a sa comprehension. D'apres 
1' invention, le modele de deplacemerit et de deformation- 
qU T on choisit pour I'objet E est affine, compose par exemple 
d'une combinaison variable en cours d f acquisition de 
translations, de rotations et d'homotheties, conformement a 
la formule (1) 



Xo = r x bcl = A x x + B x » 



a 21 (X) a 22 (X)J X + Lb 2 (^)J (1) 



ou x 0 est la position vectorielle du point P par 
rapport a une reference telle que le point O, a 1* instant de 
reference choisi pour mener les calculs et la reconstruction 

de 1' image, x est la position du point. P a un autre instant 
et notamment un instant de prise de vue, et les coefficients 
a et b dependent du temps. La partie lineaire affine du 
modele analytique du modele analytique particulier peut 
correspondre & 1 1 approximation locale affine d'une famille 
de trajectoires passant par des points situes au voisinage 
de chaque point, A. est un parametre general fdu proced£, qui 
varie entre une minimale X min et une valeur maximale 



i er depot 
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X ( event ue 1 1 ement par X lui-meme) . L 1 instant de 
reference ou 1 1 image est reconstruite correspond a A,= 0. 
La propriete d 1 attenuation ou d ! emission a calculer est 
notee f, et pour le point P qu'on a reper£ par le 

5 vecteur x on la note f x (x) a un instant X que Icon que et 

fo(xo) S. 1 1 instant de r§f6rence, o\X x est not^ xo - 

La valeur d'une projection est donnee par 

la formule (2) , ou a est un vecteur unitaire, et on 
note cette projection en fonction de ce vecteur 

10 unitaire a et du parametre X : 

g x (k, a) = Jdt f^(X) + ta ) ( O 

to* 

* II est habituel d'appliquer un filtre de 

Hilbert aux projections dans les methodes analytiques, 
ce qu'on fait ici aussi, et les projections filtr6es 
15 sont not^es selon la formule (3) 



20 



25 



- &i&(X, n ) = - J da h H (n.a)g x (X, a) , ( ^ ) 

Ou h H est donne par la formule (4) 



h H (s)=- Jdaisign(a)e 2inos . ( (j> J . 

—CO 

Les inventeurs ont etabli en outre la 
relation suivante (5) 

gjp a, A* T n) = | de ^ A ^ P oH (n, s = nJ\<a(A.))) , ("Sj 
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X max ; la position de la source S est reperee par le 
vecteur a (A,) et le temps par une fonction monotone de 
X ( e vent ue 1 1 ement par X lui-meme) . L 1 instant de 
reference ou 1 1 image est reconstruite correspond a A,=0, 
5 La propriety d 1 attenuation ou d 1 Emission a calculer est 
notee f, et pour le point .P qu'on a rep£re par le 

vecteur x on la note ^.(x) a un instant X quelconque et 

f 0 (x 0 ) a 1 1 instant de reference, ou x est note Xo . 

La valeur d'une projection est donnee par la 
10 formule (2) , oft. a est un vecteur unitaire, et on note 
cette projection en fonction de ce vecteur unitaire a 
et du parametre X : 

g x (x,a)= fdtf x (a(A)+ta) * (2) 



15 



20 



25 



II est habituel d'appliquer un filtre de 
Hilbert aux projections dans les methodes analytiques, 
ce qu'on fait ici aussi, et les projections filtrees 
sont notees selon la formule (3) 



<3*x (X, n ) =- J da b^(n.a)g x (X f a) , (3 ) 

s 1 



Oil S 1 est la sphere - unite en bidimensionnel, 
h H est donn£ par la formule (4) 

h H (s)=- Jda isign(a)e 2iTOS . (4) 



8 



I i. ou P oH designe jes projections paralleles 



V 

n-^Str-p-as — eiaiT") filtrees *par le filtre de Hilbert 
5 d'aprds la formule suivante (6) 



P 0H (n, s) = |ds' P c (n, s' )h H (s - s' ) . 



La formule d' inversion pratiquee est alors 
10 la formule suivante (7) 



f 0 (Ko) = jdX 



- r x (a(x)J 



g^(X, A x n ) 



(8) 



15 



oil est donne dans la formule suivante 

% FX 



4' A $ - 1; i det a *>hJ*' A I») } - £ < 



det A 



rn— \ 



~~ ~~ Cctte derivation s'effectue suivant le param&re X de la trajectoire du foyer, Elle est 
sp6cifiquement adaptee au modfele affine de deformation et de mobility car elle s'applique sur 
une direction Alii en considerant que la direction n est maintenue constante. 
Tl s'agit done de niaintenir constant non pas la direction orfhogonale & un rayon acquis^ mais 
Indirection oi1ho^Qnale,an.ra-v:on,6.a uivalent de lajg ^om^n^g^ve— : — — 

— ^3jnestla direction orthogonale k la droite de projection & Finstant X passant par le point 

mobile correspondant au point P. 
1 A Cstte formuie d ! inversion tient compte des 

deformations de l'objet E et comporte, par rapport a 

d'autres formules, etablies dans des cas ordinaires, 

une pond6ration en fonction de la deformation de 

l'objet (par le determinant de A x ) et de la position de 

25 la trajectoire (par la norme entre x 0 et T x (a(?i)). Ces 
conditions peu contraignantes s'expliquent parce que, 
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Les inventeurs ont etabli en outre la relation 
suivante (5) 



|detA x | 



P oH (n,s = nT x (a(X)), 



(5) 



10 



ou P oH designe des projections paralleles P Q 
(fi,s), obtenues par rearrangement des projections 
ef fectivement mesurees sur l'objet a 1 ' etat de 
reference et filtrSes par le filtre de Hilbert d'apres 
la formule suivante (6) 



PaH(n,s)= ^ds'P^s^hHCs-s') . 



15 



La formule d* inversion pratiquee est alors la 
formule suivante (7) 



J 1 



T — * 

•g F x(X,A A n ) 



'.A r (*b) |^o-r x (a(A)| 
oil gpx est donne dans la formule suivante (8) 



g^^Ajn)=|-jdet A x |g HX (x, Ajn) }- JL 



detA. 



Ain 



•9HX 



Ain 



2 0 Cette derivation s'effectue suivant le 

parametre X de la trajectoire du foyer, Elle est 
specif iquement adaptee au module affine de deformation 
et de mobilite car elle s» applique sur une direction 
A*fi en considerant que la direction n est maintenue 

25 constante. II s 1 agit done de maintenir- Constant, non 
pas la direction orthogonale k un rayon acquis, mais la 



i er uepoi 



ivioainee le ^3/ 



avec un espace affine de transformation, les droites de 
projection de la figure 1 restent des droites dans la 
geometrie fictive de la figure 2, de sorte que le 

, iwj: J*??* . 

probleme numerique r^ I ^s^pars^e^^eeup— GeH^l-a-gue^ 7 
Dans la formule (7) 

bornes de 1 1 integral e A r (x 0 j designent un ensemble 

minimal de positions X de la source S tel que sur 
l'objet E a 1 1 etat de reference les directions des 

droites reliant xo k r x (a(k)) couvrent tout l'intervalle 
d'un demi-tour de trajectoire sans redondance . Si on 
prend soin de limiter ainsi l'int^grale, la formule de 
reconstruction est parfaite, 

II est manifeste que les deformations 
af fines ne suffisent pas a mod61iser tous les objets 
qu'on peut devoir gtudier dans la pratique. Refaisant 
le raisonnement des figures 1 et 2, on donne les 
figures 3a , 3b , 3c, et 4. L'objet E evolue de fagon 
complexe entre les positions indiquees successivement 
aux etats Xl, \2 et %3 oil la position de la source S a 
aussi #te indiquee de mime que le rayon Rl, R2 ou R3 
menant au point P. 

En remettant l'objet E dans un etat de 
reference, la figure 4 montre que les rayons notes R'l, 
R f 2 et R ! 3 menant au point P et correspondant aux 
rayons Rl, R2 et R3 (passant par les mSmes points de 
l f objet E) ne sont plus rectilignes. La formule (7) 
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direction orthogonale au rayon equivalent de la 
g^ometrie fictive. 

Ala est la direction orthogonale a la droite de 

projection & 1 1 instant X passant par le point mobile 
5 correspondant au point I?. 

Cette formule d 1 inversion tient corapte des 
deformations de I'objet E et comporte, par rapport a 
d 1 autr.es formules, etablies dans des cas ordinaires, 
une ponderation en fonction de la deformation de 
10 I'objet (par le determinant de A x ) et de la position de 

la trajectoire (par la norme entre x 0 et r x (a(X)). Ces 
conditions peu contraignantes s'expliquent parce que, 
avec un espace affine de transformation, les droites de 
projection de la figure 1 restent des droites dans la 
15 .geometrie fictive de la figure 2, de sorte que le 
probldme numerique peut §tre rgsolu analytiquement . 
Dans la formule (7) 



bornes de 1 1 integrale A r (x 0 ) d^signent un ensemble 
20 minimal de positions X de la source S tel que sur 
I'objet E a 1 1 <§tat de reference les directions des 
droites reliant x 0 & I\(a(A,)) couvrent tout 1 'interval! e 
d'un demi-tour de trajectoire sans redondance . Si on 
prend soin de limiter ainsi 1 ' integrale, la formule de 
25 reconstruction est parfaite. 

II est manifeste que les deformations affines 
ne suffisent pas a modeliser tous les objets qu'on peut 
devoir etudier dans la pratique. Refaisant le 
raisonnement des figures 1 et 2, on donne les figures 



f 0 (x 0 ) * J dX. pr: z— ug FA (^ A x n ) n'est alors 

af ||x 0 - r x (a(x)\\ 

plus valable directement. On s*y ramSne toutefois par 
des approximations , en supposant une deformation affine 
de l'objet E particuliSre a chaque point P et valable 
5 autour de lui* On calcule alors la matrice de 
deformation correspondant & A x de la formule (1) 



10 



15 



20 



par la formule (9) 



A x x + 



_ ra lx (A.) a 12 (*.)] - [-^(X)! 
" kx(A.) a 22 (*)J X + |ba(*)J 



- r Xi (X, r _1 (x 0 ) -^pj^f$ 



dx 



25 



oil r xl (x f xjet r x2 (A,, x) sont les composantes 

suivant les directions xl et x 2 principales de r(x,, x) . 
BX. n'a pas besoin d'etre calcule. les deplacements 
estimes de tous les points de l'objet E sont done les 
coefficients r(x,, x) qui peuvent concrdtement Stre 
introduits dans des cartes de deplacement qui sont lues 
et exploitees au moment de I 1 inversion. 

A chaque point P, il faudra encore 

determiner 1 1 ensemble de projections A r (xo) qui seront 
utilis^es pour l 1 inversion, selon le principe precedent 
qu'un intervalle angulaire d f un demi-tour doit etre 
couvert par ces projections. Comme les rayons de la 
ggometrie fictive R'l, R'2 et R'3, etc. ne sont plus 
rectilignes, on considere plutdt les directions de leur 
tangente au point P d 1 intersection. 
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raisonnement des figures 1 et 2, on donne les figures 
3a, 3b, 3c, et 4. L'objet E evolue de fagon complexe 
entre les positions indiquees successivement aux etats 
Xl, X2 et Xz oH la position de la source S a aussi £te 
indiquee de meme que le rayon Rl, R2 ou R3 tnenant au 
point P. 

En . remettant l'objet E dans un Stat de 
reference, la figure 4 montre que les rayons notSs R'l, 
R. ! 2 et R f 3 menant au point P et correspondant aux 
rayons Rl, R2 et R3 (passant par les m§mes points de 
l'objet E) ne sont plus rectilignes. La formule (7) 



Co(xo)- JdA 
Ar 



p-r x p(x)| 



g FX (A,A£n ) n'est alors plus 



valable directement. On s'y ramene toutefois par des 
approximations, en supposant une deformation affine de 
l f objet E particulidre a chaque point P et valable' 
autour de lui. On calcule alors la matrice de 
deformation correspondant a . de la formule (1) 



*o = T x (x) - A^x + B x = 
la formule (9) 



a 21 (A,) a 22 (X) 



)J b>M} 



par 



oil r^A,, x) et T^ffa, x) sont les composantes 

suivant les directions xl et x 2 principalis de T^k, x) • 
BA, n ! a pas besoin d'Stre calcule. Les deplacements 



tvioaiTiee le zj/a 
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La formule d 1 inversion peut alors s ' ecrire 
selon la formule (10) 



/. J$) « f dX -j- l—-^^ t n(T-^ 0 ) X) ) , 

0 ATg.) ii xo - r x fa(xjl 



oil 



5 ggfc est exprime par la formule (11) 

det A x (r~ x (x 0j 



9sx (x, n(r-^ 0 } x)) = ± 



►g^X, n^-^xo) x)) 



10 



15 



De plus, n(r _1 (xo) x) est la direction 
orthogonale a la droite acguise a 1 ' instant X passant 
par r~(x 0 ), c'est-a-dire la retroprojection de la 
projection passant par le point ; ^(t^xo) x)) est un 
facteur lie & la deformation des longueurs sur les 
rayons, d'apres la formule (12) 




Enfin, I^xo - r x (a(XjJ est la distance du point P a 1 « image 

de la source sur l'objet de reference, le long du rayon 
20 virtuel R' . 

On s'est interess<§ jusqu'ici a des 
reconstructions de l'objet E par tranches, dans des 
conditions en eventail parallele. Le procede precedent 
peut etre etendu avec des acquisitions par un 
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10 



15 



estim£s de tous les points de I'objet E sont done les 
coefficients r{x, x) qui peuvent concretement §tre 
introduits dans des cartes de deplacement qui sont lues 
et exploit£es au moment de l 1 inversion. 

A chaque point P, il faudra encore determiner 

1' ensemble de projections A r (xo) qui seront utilises 
pour 1' inversion, selon le principe precedent qu'un 
intervalle angulaire d'un demi-tour doit §tre couvert 
par ces projections. Comme les rayons de la g^ometrie 
fictive R'l, R'2 et R ! 3, etc. ne sont plus rectilignes, 
on considdre plutot les directions de leur tangente au 
point P d 1 intersection. 

La formule d 1 inversion peut alors s'ecrire 
selon la formule (10) 



-4 



est exprime par la formule (11) 



^5> x (A,n(r 1 ^o)x)) , ou g P x 



20 



25 



det AjJ 


m 







V^X, n(r^\x)) 



De plus, n^^xojlx) est la direction orthogonale 

a la droite acquise a 1 1 instant X passant par r^xo) / 
e'est-a-dire la retroprojection de la projection 
passant par le point ; pfn^^xo), x)) est un facteur li<§ & 
la deformation des longueurs sur les rayons, d'apres la 
formule (12) 
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rayonnement conique dans les conditions class iques de 
la figure 5, ou yobjet E est represents dans son 
entierete tridimengionnelle; de meme que le reseau de 
detecteurs D, qui est compose d'une serie de couches de 
detecteurs semblables a celle des figures precedentes ; 
un rayon R peut etre exprime par trois parametres, a 
savoir celui X de la position de la source S, un angle 
Y que fait le rayon R dans le plan de la trajectoire T 
par rapport a 1 • axe central X du faisceau, et une cote 
q reperant la couche de detecteurs a laquelle le rayon 
aboutit. L 1 angle y remplace le vecteur de direction S 
precedemment utilise dans les formules par commodite. 
La formule suivante (13) 

\ •£ (*.))) 

donne 1 • inversion qu'on entreprend alors 
dans ces conditions geometriques particulieres et en 
faisant une approximation appelee de Feldkamp courante 
dans l'art pour traiter les projections coniques . Les 
mesures sont ici multipliees par cos A, qui est la 
ponderation classique compensant le caractere oblique 
du rayon acquis. Cet angle A est donne sur la figure 5 
et represente 1' angle du rayon R avec le plan de la 
trajectoire T. De plus, l(x 0 - r x (a(xj est la distance de 
x a 1' image de la source S sur l'objet E a 1 ' etat de 
reference le long du rayon virtuel R- , en projection 
sur le plan de la trajectoire T, et n (^(xo) x) est la 
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Enf in, l(x 0 - r x (a(k)) est la distance du point P a 
l 1 image de la source sur I'objet de reference, le long 
du rayon virtuel R f . 
5 On s'est interessS jusqu'ici & des 

reconstructions de l'objet E par tranches, dans des 
conditions en <§ventail parallele. Le procSde pr<§c§dent 
peut §tre Stendu avec des acquisitions par un 
rayonnement conique dans les conditions class iques de 

10 la figure 5, ou I'objet E est represents dans son 
entierete tridimensionnelle, de meme que le reseau de 
detecteurs D, qui est compost d'une serie de couches de 
detecteurs semblables & celle des figures pr6c£dentes ; 
un rayon R peut etre exprime par trois parametres, a 

15 savoir celui X de la position de la source S, un angle 
y que fait le rayon R dans le plan de la trajectoire T 
par rapport & I'axe central X du faisceau, et une cote 
q reperant la couche de detecteurs a laquelle le rayon 
aboutit. L 1 angle y remplace le vecteur de direction a 

20 pr^cedemment utilise dans les formules par commodite. 
La formule suivante (13) 

25 donne 1 1 inversion qu'on entreprend alors dans ces 
conditions geometriques particulieres et en faisant une 
approximation appelSe de Feldkamp courante dans 1 1 art 
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direction orthogonale a la droite acquise a 1' instant X 
passant par T~ (x 0 ) sur ce plan. Enf in, on a vu que 
q(k, r" 1 ^)) est la cote sur le reseau de dStecteurs D de 

la droite,, acquise a 1 1 instant X, passant par' x 0 ♦ 
5 Un procede particulier, dScrit dans le 

brevet frangais d'abord d#pos§ sous le numero 01 07918, 
peut encore §tre applique ici moyennant quelques 
transformations. Ce procede s T applique a des 
acquisitions rapides menses autour de l'objet en 

10 faisant plusieurs tours et consiste en une 
reconstruction de l'objet -e^t^^^^r#el en evaluant son 
contenu par des blocs de projections prises sur tine 
fraction de tour seuleraent de la source. Les inversions 
de projections menees sur les blocs produisent des 

15 sous-images qui sont incorrectes puisqu'elles ne 
comprennent qu'une partie des raesures, mais sont 
obtenues tres rapidement et permettent d'obtenir avec 
une m§me rapidite la loi de deplacement ou de 
deformation de l'objet en comparant des sous images 

2 0 homologues, prises pour des mgraes positions de sources 
espac€es d'un tour complet ou d'un demi-tour ; les sous 
images sont reconstruites a un instant de reference 
pour chacun des groupes de blocs et finalement 
combinees entre elles pour dormer 1 r image complete de 

25 l'objet. Ici, les limites des blocs sont donnees dans 
la g^ometrie virtuelle. Le procede n'est par ailleurs 
modif ie. ^ ^ 

Enf in, il est possible d'ameliorer le 
procede pour 1' etude d'un objet soumis a des phenomenes 
30 periodiques. On select ionne alors comme blocs d' etude 
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pour traiter les projections coniques. Les mesures sont 
ici multipliees par cos A, qui est la ponderation 
classique compensant le caractSre oblique du rayon 
acquis. Cet angle A est donne sur la figure 5 et 
5 represente l 1 angle du rayon R avec le plan de la 

trajectoire T. De plus, Lp 0 - r x (a(jC)) est la distance de 

x a 1 1 image de la source S sur l'objet E a I'etat de 
reference le long du rayon virtuel R 1 , en projection 

sur le plan de la trajectoire T, et n (r'^xo), x) est la 
10 direction orthogonale a la droite acquise a 1 1 instant X 
passant par T^Cxo) sur ce. plan. Enfin, on a vu que 
q(X, r z -1 (x 0 )) est la cote sur le reseau de detecteurs D de 
la droite, acquise & l 1 instant X, passant par le point 
mobile correspondant k x 0 . 
15 Un procede particulier, d^crit dans le brevet 

fran<?ais d'abord depos6 sous le numero 01 07918, peut 
encore Stre applique ici moyennant quel que s 
transformations. Ce procede s 1 applique a des 
acquisitions rapides menees autour de l'objet en 

2 0 faisant plusieurs tours et consiste en une 

reconstruction de l'objet anime en £valuant son contenu 
par des blocs de projections prises sur une fraction de 
tour seulement de la source. Les inversions de 
projections menees sur les blocs produisent des sous- 
25 images qui sont incorrectes puisqu'elles ne comprennent 
qu'une partie des mesures, mais sont obtenues trSs 
rapidement et permettent d'obtenir avec une meme 
rapidite la loi de deplacement ou de deformation de 
1' ob jet en comparant des sous images homologues, prises 

3 0 pour des m§mes positions de sources espacees d f un tour 
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seulement ceux qui sont pris a des phases du phenomene 
analogues a celles de 1' instant de r6f£rence. 
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complet ou d'un demi-tour ; les sous images sont 
reconstruites Sl un instant . de reference pour chacun des 
groupes de blocs et f inalement combinees entre elles 
pour dormer 1 1 image complete de 1 r ob j et . Ici , les 
5 limites des blocs sont donnees dans la geometrie 
virtuelle. Le proced£ n'est pas par ailleurs modifiS. 

Enfin, il est possible d'ameliorer le precede 
pour l ! £tude d'un objet soumis a des phenomdnes 
pSriodiques. On selectionne alors comme blocs d'^tude 

10 seulement ceux qui sont pris & des phases du phenomene 
analogues & celles de I 1 instant de reference. 

On pourra se reporter aux organigratnmes des 
figures 6, 7 et 8 pour finir de prendre connaissance de 
1" invention. lis detaillent trois modes de realisation 

15 de l 1 invention, ou la reconstruction est faite 
respect ivement avec compensation du mouvement ; 
compensation du mouvement et compensation temporelle ; 
et compensation du mouvement, compensation temporelle 
et prise en compte d'une pSriodicite de Involution de 

20 1» objet. Ce qui a ete decrit auparavant est applique 
dans ces procedes . Les voxels* sont evidemment les 
points de 1 1 image de 1 1 ob j et qu 1 on considSre dans 
1 1 image reconstruite . 
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complet ou d'un demi-tour ; les sous images sont 
reconstruites a un instant de reference pour chacun des 
groupes de blocs et finalement combinees entre elles 
5 pour donner 1' image complete de l 1 objet. Ici 7 les 
limites des blocs sont donnees dans la geometrie 
virtuelle. Le procede n'est pas par ailleurs modifie. 

Enfin, il est possible d'ameliorer le procede 
pour 1' etude d'un objet soumis a des phenomenes 

10 periodiques. On selectionne alors . corame blocs d T etude 
seulement ceux qui sont pris & des phases du phenomene 
analogues a celles de 1' instant de reference. 

Si le rayonnement est conique, le procede peut 
comprendre, avant l'etape de filtrage des mesures 

15 acquises, une etape de ponderation des mesures acquises 
adaptees a une acquisition des mesures en geometrie 
conique ; la retroproj ection est alors effectuee selon 
la geometrie conique en tenant compte du modele 
analytique de mobilite et de deformation de 1' objet. 

20 On pourra se reporter aux organigrammes des 

figures 6, 7 et 8 pour finir de prendre connaissance de 
1' invention. lis detaillent trois modes de realisation 
de 1 T invention, ou la reconstruction est faite 
respectivement avec compensation du mouvement ; 

25 compensation du mouvement et compensation temporelle ; 
et compensation du mouvement, compensation temporelle 
et prise en compte d'une periodicite de 1« evolution de 
1 ! objet. Ce qui a ete decrit auparavant est applique 
dans ces procedes. Les voxels sont evidemment les 

30 points de 1 1 image de 1' objet qu'on considere dans 
1 1 image reconstruite . 
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RE VEND I CAT I ONS 

1) Procede de reconstruction d'une image 
toraographique d f un objet notamment mobile et 
deformable, 1 1 image etant un ensemble de valeurs d'une ^*$C<; 
propriete prises par 1' objet a une collection de points 1/**** 
de l 1 objet, comprenant l'emploi : d'un rayonnement P*. 
divergent & partir d'un foyer et traversant I 1 objet, le / s 
foyer etant mobile autour de 1 ! objet ; d'un modele^ <* 
analytique de mgbilite et de deformation de 1' objet ; 
et d'unfc #©aanule analytique pour obtenir lesdites jJr aj^ 



valeurs S partir de sommes des valeurs de la propriete 
le long de lignes de projection menant au foyer et 
passant respectivement par les points ; caractSrisS eny J^*^ 
ce que le mobile est affine, comprenant des 4 "**^v^* V 
translations, des rotations et des homotheties de... % 0y j^ ^ 
1' objet £ partir d'une origine, et en ce que 1^43XDEUiJe..v 

^ = Itu^rM^ Al ' • 



(?) 



2) Procede de reconstruction d'une image 
tomograph! que selon la revendication 1, caracterise en 
ce qu'il comprend une etape de definition d'intervalle 

d 1 integration donnant f 0 (xo) de fagon que les lignes de 

projection aient des directions s'etendant sur un 
it K) 

intervalle d^xi demi-tour5 considers dans un Stat de 
1' objet ou 1 1 image est reconstruite 

3) Procede de reconstruction d'une image 
selon I'une quelconque des revendi cat ions 1 ou 2, 
caracterisS en ce qu'il est applique avec un modele 
analytique de mobilite et de deformation particulier §l 
chaque point de 1' objet. 
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REVEND I CATI ONS 

1. Precede de reconstruction d'une image 
tomographique d'un objet notamment mobile et 
5 deformable, 1 1 image etant un ensemble de valeurs d'une 
propri£te prises par 1 'objet £ une collection de points 
de 1» objet, comprenant l'emploi : d ! un rayonnement 
divergent a partir d'un foyer et traversant l 1 objet, le 
foyer 6tant mobile autour de l l objet ; d'un moddle 

10 analytique de mobilite et de deformation de 1' objet 
defini pour chaque position du foyer ; et d'un procede 
de calcul analytique pour obtenir lesdites valeurs & 
partir de sommes des valeurs de la propriete le long de 
lignes de projection menant au foyer et passant 

15 respect ivement par les points ; caracterise en ce que 
le mobile est affine, et est une combinaison variable 
en cours d 1 acquisition, cette combinaison, comprenant 
des translations, des rotations et des homotheties de 
1 'objet a partir d'une origine, et en ce que le procede 

20 . de calcul analytique comprend les etapes suivantes : 

- f iltrage des mesures acquises par un f iltre de 
Hilbert ; 

- pond^ration des mesures filtrees, cette 
pond§ration etant d6pendante du moddle 

25 analytique de mobilite et de deformation de 

1 1 ob j et ; 

- derivation des mesures filtrees ponderees 
suivant la trajectoire du foyer en considerant 
une direction adaptee au modele, cette 

3 0 direction etant maintenue constante, et 

obtention de mesures modifi^es ; 



25 



2^ 

16 ^W""* *^ 

4) ProcSde^ de reconstruction d'une image ^ 
selon la revendication 3, caract^rise en ce que _3re 
module analytique particulier est: SaTcul£ par la— 



5) Proced^ de reconstruction d'une image 
selon les revendications 2 et 3, caracterisS en ce que 
les direction'des lignes de projection sont considerees 
au point de l'objet oil le proced£ est appliqu6. 

10 6) Procede de reconstruction d'une image 

selon 1'une quelconque des revendications 1 §. 5, 
caract^risS en ce que la reconstruction s'effectue 
d'abord par sous-images construites avec des lignes de . 
projections ayant des directions comprises dans les. .' 

15 intervalles angulaires plus petits gu'un demi-tour, etr- - 
en ce qu'il comprend une 6tape de determination des * 
intervalles angulaires auxquels appartiennent les 

lignes de projections .-^-> v~ c*. ^ ^Ifl c^j; n \ 

g JJf Procede de reconstruction d'une image / 

2 0 selon l'une quelconque des revendications precedentes r 
ou le rayonnement est conique, earracrtrgr drser en que !ba- 



JA _*.. r r . . 



— - : : — 7 caracttrise en ce que la loi devolution g&i&ale de Pobjet 

est penodique et en ce que les instants de r6fSrence sont choisis pour la meme phase des 
JB*!SS^ 1'pbjet (Cf^reiign^ 
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- r^troprojection des mesures tnodifiees. 

2. Procede de reconstruction d'une image 
tomographique selon la revendication 1, caractSrisli en 

5 ce qu'il comprend une etape de definition d'intervalle 

d 1 integration dormant fo(x 0 ) ^ e fagon que les lignes de 
projection aient des directions s 1 Stendant sur un 
intervalle de N demi- tours consider^ dans un 6tat de 
l'objet ou 1 T image est reconstruite, N gtant sup^rieur 
10 ou egal a 1. 

3. ProcedS de reconstruction d'une image selon 
l'une quelconque des revendications 1 ou 2, caracteris^ 
en ce qu'il est applique avec un modele analytique de 

15 raobilite et de deformation particulier 3. chaque point 
de l'objet. 

4. Procede de reconstruction d'une image selon 
la revendication 3, caracterisS en ce que la partie 

20 lineaire affine du modele analytique particulier 
correspond a 1 1 approximation locale affine d'une 
famille de trajectoires passant par des points situes 
au voisinage de chaque point, 

25 5. Procede de reconstruction d'une image selon 

les revendications 2 et 3, caract^rise en ce que les 
directions des lignes de projection sont considerees au 
point de l'objet ou le procedS est applique. 

3 0 6. Procede de reconstruction d'une image selon 

l'une quelconque des revendications 1 a 5 f caracteris§ 
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- filtrage des mesures acquises par up filtrc de Hilbert } 

- pond&ration des mesures filtrSes , cette pondfoation etant dependante du module 
anaiytique de mobilite et de deformation de T objet {€tyS ^6s 2^^ j 

- derivation des mesures filtr6es ponderees sitivant la trajectoire du foyer en considcrant 
une direction adaptSe au module (Cfr^ligae-l 7~a-2(Hnedifiee£), cette direction Slant 
maintenuc constante et obtention de mesures modifiees * 

- r&roprojection des mesures modifies. ) 



ftev-6^i*4ti^om^ quia^omprendrune etape de 

combinaison des sous-images apr&s carnation de la loi d'evolution generate de I'objet 
jtesc^iption p. 3-ligne 6-3-2) entre les positions associees a diff&rents instants de reference 
assoctes & chaque intervalle angulaire. £ ^ ^ } 



caracterise en ce qu'avant T6tape de filtrage des mesures acquises, le procede comprend une" 
&ape de pond6ration des mesures acquises adarptee a unc acquisition des mesures en 
g6om&rie conique et en ce que la r6troprojection est effectu<Se selon. la geom^tricconique en 
tenant compte du module anaiytique de mobilite et de deformation de Pobjet # 4* j 
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en ce que la reconstruction s'effectue d'abord par 
sous-images construites avec des lignes de projections 
ayant des directions comprises dans les interval les 
angulaires plus petits qu'un demi-tour, en ce qu'il 
comprend une etape de determination des intervalles 
angulaires auxquels appartiennent les lignes de 
projections et en ce qu*il comprend une £tape de 
combinaison des sous- images aprds correction de la loi 
devolution g^n^rale de I'objet entre les positions 
associees a differents instants de reference associes & 
chaque intervalle angulaire. 

7. Proc£d£ de reconstruction d'une image selon 
i'une quelconque des revendi cat ions prec^dentes , 
caracterisS en ce que la loi devolution generale de 
I'objet est pgriodique et en ce que les instants de 
r£fSrence sont choisis pour la m§me phase des periodes 
de mouvement de 1 1 ob j et . 

8, ProcedS de reconstruction d'une image selon 
l'une quelconque des revendi cat ions precedentes, oil le 
rayonnement est conique, caracterise en ce qu'avant 
l'§tape de filtrage des mesures acquises, le procede 
comprend une etape de pondgration des mesures acquises 
adapt§e a une acquisition des mesures en ggom^trie 
conique et en ce que la rgtroprojection est effectuee 
selon la geometrie conique en tenant compte du module 
analytique de mobilite et de deformation de I'objet. 



• 1/// • 
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FIG. 5 






Acquisition des projections 
d'un objet evoluant sur une 
duree d&erminee 



1A 



synchronisation 



IB 



Choix d'un instant 
de reconstruction 



ID 



Acquisition des trajectoires de 
tous les voxels sur la duree 
d' acquisition (C£ &T 



1C 



Calcul des approximations locales 
affines en tout voxel a tout instant par 
rapport a l'instant de reconstruction 



T 



Selection des ensembles de mesures 
optimaux pour la reconstruction 
specifique a chaque voxel et 
dependant de sa deformation locale 



IE 



IF 



Reconstruction adaptee a la d6formation 1 
h un instant donnS de chaque voxel de 
l'objet 



I 



Image reconstruite avec 
reduction des artefacts 



1G 



1H 



Auteur :..S(£ bastien Roux 
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Acquisition des projections d'un 
Dbjet 6voluant sur une duree 
determinee 



synchronisation 



4B 




Choix d'un instant 
de reconstruction 



4D 



Acquisition des trajectoires de 
tous les voxels sur la dur6e 
d' acquisition 



4C 



Calcul des approximations locales 
affines en tout voxel a tout instant par 
rapport k 1'instant de reconstruction 



4E 



Selection des groupes de blocs de 
mesures optimaux pour la 
reconstruction partielle sp6cifiques a 
chaque voxel et dependant de sas 
deformation locale 



4F 



Reconstruction adapt6e & la d6formation 
k un instant donn<5 de chaque voxel de 
l'objet et avec compensation des 
variations temporelles (C££$ 



I 



4G 



Image reconstitute avec 
r6duction des artefact? 



4H 



7 
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Acquisition des projections d'un 
objet evoluant sur une duree 
determinee 



synchronisation 



6B 



Choix d'un instant 
de reconstruction 



6D 



Signal synchronisateur avec 
la phase de revolution 
p&iodique 



Acquisition des trajectoires de 
tous les voxels sur la duree 
d' acquisition (CWf 



6F 



6C 



Calcul des approximations locales 
affines en tout voxel h tout instant 



6E 



S61ection des groupe de blocs 
ensembles de mesures optimaux de 
meme phase que Pinstant de 
reconstruction pour la 

reconstruction partielle sp6cifique h. 
chaque voxel et dependant de sa 
deformation locale : 

~ 6G 



Reconstruction adapt^e a la deformation 
a un instant donne de chaque voxel de 
Tobjet et avec compensation des 
variations temporelles (Cf 5) 



6H 

▼ 

Image reconstruite avec 
reduction des artefact. 
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1 S'H v a olus de trois inventeurs, utilisez plusieurs formulaires. Indiquez en haut a droite le N° de la page sum du nombre de pages. 


I DATE ET SIGNATURES) 

DU (DES) DEMAND EUR(S) 
1 OU DU MANDATAIRE 
1 (Nom et quality du signataire) 

1 PARIS LE 23 JU1LLET 2003 
I J. LEHU 





Elle garantit un droit d'acc&s et de rectification pour les donnees vous concernant annrfcc Ha imnpi 




DfcPARTEMENT des brevets 

26 bis, rue de Saint Petersbourg 
75800 Paris Cedex 08 

Telephone : 33 (1) 53 04 53 04 Telecopie : 33 (1) 42 94 86 54 



regue le 23/07/03 



BREVET DWVENT10N 
CERTIFtCAT D'UTILITE 

Code de la propriete intellectueUe - Uvre VI 

DESIGNATION D'INVENTEUR(S) Page N* ?. . /? 

(A fournir dans le cas 0C1 les demandeurs et 

les inventeurs ne sont pas tes m§mes personnes) 

Cet ImprlmS est & remptir llsiblement a 1'encre noire 



« 



N° 11235 



DB113@W / 27 



Vos references pour ce dossier (facultatif) 



B 14364.3/JCI 



N° D'ENREGISTREMENT NATIONAL 



03.07848 DU 27.06.2003 



T1TRE DE L'INVENTION (200 caracteres on espaces maximum) 

PROCEDE DE RECONSTRUCTION D'UNE IMAGE TOMOGRAPHIQUE PAR UNE METHODE ANALYTIQUE 
COMPRENANT UNE MODELISATION AMELIOREE DU MOUVEMENT DE L'OBJET. 



LE(S) DEMANDEUR(S) : 

COMMISSARIAT A L'ENERGIE ATOMIQUE 
31-33 rue de la Federation 
75752 PARIS 15 erne. 

CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE 
3 rue Michel Ange 
75794 PARIS CEDEX 16 

DESIGNE(NT) EN TANT QU'INVENTEUR(S) : 



Q Nom 


GRANGEAT 


Prenoms 


Pierre 


Adresse 


Rue 


133, Altee Vignes du persan 


Code postal et ville 


(3 1 8 13 13 10 I SAINT ISMIER 


Societe d'appartenance (facultaiij) 




Q Nom 




Prenoms 




Adresse 


Rue 




Code postal et ville 


I I I I I I 


Societe d'appartenance (facultatif) 




O Nom 




Prenoms 




Adresse 


Rue 




Code postal et vide 




Societe d'appartenance (facultatif) 





S'il y a plus de trois inventeurs. utilisez plusieurs formulaires. Indique2 en haut a droite le N° de la page sutvi du nombre de pages 



DATE ET SIGNATURE(S) 
DU (DES) DEMAND EUR(S) 
OU DU SWAN DATA! RE 
(Nom et qualite du signataire) 

PARIS LE 23 JUILLET 2003 
J. LEHU 




La loi n°78-17 du 6 janvier 1978 relative a rinformatique, aux fichlers et aux liberies s'applique aux reponses faltes a ce fbrmulaire. 
Elle garantit un droit d'acces et de rectification Dour les donnees vous concsrnant aimr^ Ha imnpi 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 
^BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

"^hLINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



